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ABSTRACT H 

Wave rotors, used in a gas turbine topping cycle, offer a i 
potential route to higher specific power and lower specific fuel 
consumption. In order to exploit this potential properly, it is L 

necessary to have some realistic means of calculating wave M, 

rotor performance, taking losses into account, so that wave n 

rotors can be designed for good performance. This in turn P, 

requires a knowledge of the loss mechanisms. The experiment Pr ( 

reported here was designed as a statistical experiment to Tj 

identify the losses due to finite passage opening time, friction, Tr, 

and leakage. For simplicity, the experiment used a 3-port, flow t 

divider, wave cycle, but the results should be applicable to U 

other cycles. A 12" diameter rotor was used, with two different 
lengths, 9" and 18", and two different passage widths, 0.25" V, 

and 0.54”, in order to vary friction and opening time. To vary Wj 

leakage, moveable end-walls were provided so that the rotor to a 

end-wall gap could be adjusted. The experiment is described, (3 

and the results are presented, together with a parametric fit to 
the data. The fit shows that there will be an optimum passage y 
width for a given wave rotor, since, as the passage width 6 
increases, friction losses decrease, but opening-time losses 5 P 

increase, and vice-versa. Leakage losses can be made small at 
reasonable gap sizes. 5 2 (x) 

s 


NOMENCLATURE 


a 

B 

bo,b„bii 


C P 

D h 

Z>(x) 

F 

G, G p , G„ 


speed of sound 

width of a passage on the rotor 

constants defined in Eq. (1 9) 

specific heat at constant pressure 

hydraulic diameter of passages 

drag on a plate of length x 

dimensionless friction parameter, defined in Eq. (4) 

leakage parameters, defined in Eqs. (15), to (17) 


•'p 

0 . 

win*, 


0, 

e„ 


height of a passage on the rotor 
angle of incidence of entering air in the rotor refer- 
ence frame 
length of the rotor 
Mach number in region j 
potential number of cycles on the rotor 
absolute stagnation pressure in region j 
relative stagnation pressure in region j 
absolute stagnation temperature in region j 
relative stagnation temperature in region j 
time 

circumferential velocity of the rotor at the average 

radius of the passages 

absolute velocity in region j 

relative velocity in region j 

angle of flow in inlet duct 

ratio of mass flow in high pressure port to total 

mass flow 

ratio of specific heats 
end-wall to rotor gap spacing 
end-wall to rotor gap spacing at a port 
end-wall to rotor gap spacing away from a port 
boundary layer momentum thickness 
expansion ratio - ratio of the pressure in the low 
pressure port to the pressure in a passage just 
before reaching the low pressure port 
efficiency 

angular extent of low pressure port 
angular extent of low pressure region at inlet 
endwall 

angular extent of low pressure region at exit 
endwall 
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V 

kinematic viscosity 

p 

gas density 

X 

dimensionless opening time, defined in Eq. (1) 

0) 

rate of rotation of the rotor 

Subscripts 

cav 

cavity surrounding the rotor into and from which 
leakage occurs 

i 

general subscript for any of the three subscripts 
below 

in 

inlet port 

hi 

high pressure port 

lo 

low pressure port 

pas 

passage immediately before opening to the low 
pressure port 


INTRODUCTION 

The performance of gas turbine engines can be improved if 
the combustion temperature can be raised. This is difficult to do 
because the turbine inlet temperature is limited by material con- 
siderations (Peacock and Sadler, 1992). Increased perform- 
ance can be achieved if the pressure entering the turbine can 
be increased, while maintaining the rate of heat addition. Thus 
the combustion step should be configured so as to result in a 
pressure gain, rather than a pressure loss. Two techniques for 
achieving this are unsteady combustion (Kentfield, 1995), and 
use of a wave rotor topping cycle (Meyer, 1947, Zauner et al., 
1993, Kentfield, 1995). Because unsteady combustion currently 
shows only modest pressure gains, the wave rotor approach 
seems preferable. Calculations show that increases of 20% in 
specific power, and reductions in specific fuel consumption of 
18% are possible by using a wave rotor topping cycle (Wilson 
and Paxson, 1995). 

Wave rotors are devices which use unsteady waves to 
produce steady streams of gas which are at either higher or 
lower stagnation pressure than the input stream. The rotor itself 
has a set of passages on its periphery. In the present exper- 
iment, which was designed for no exchange of shaft work, the 
passages are straight, and aligned axially. As the rotor rotates, 
these passages are alternately exposed to ports at differing 
pressures. Typically, at the exhaust, or low pressure, port, the 
passage contains gas at some higher pressure just before the 
passage rotates into juxtaposition with the port. Exposure to the 
low port pressure causes an expansion wave to propagate into 
the passage. Later in the cycle, the passage, now at lower 
pressure, will be opened to the inlet port, where the gas is at 
higher pressure, thereby causing a shock wave to be pro- 
pagated into the passage, increasing the stagnation pressure of 
the gas. The exact sequence of waves will depend on the cycle 
employed. Several different cycles are possible, each serving a 
different function. Examples are three-port cycles, used as flow 
dividers or equalisers (Kentfield, 1963), four-port cycles, used 
for superchargers (Jenny and Zumstein, 1982), topping cycles 
for gas turbine engines (Meyer, 1947, Zauner et al., 1993), and 
a wave superheater wind-tunnel (Weatherston et al., 1959), and 
five, and nine port cycles, again intended for use as topping 


cycles (Thayer et al., 1981). In addition, wave engines for gen- 
erating shaft work have been developed (Pearson, 1985, 
Weber, 1995). However, all the cycles have common features, 
and a study of one should produce results applicable to all 
cycles, at least in a general sense. For example, all cycles 
employ an expansion wave, and expanding through too large a 
pressure ratio leads to losses for any cycle. Obviously, for 
maximum output, whatever the application, the efficiency of the 
wave rotor should be as high as possible, i.e. the losses should 
be minimised. In order to do this, it is necessary to know the 
source of the losses, and their dependence on controlling 
parameters. This study is aimed at assessing experimentally 
the magnitude of various wave rotor losses as a function of the 
parameters which affect them. This will be achieved by measur- 
ing the performance of a wave rotor as various geometrical 
parameters (passage width, rotor length, and rotor-casing clear- 
ance) are varied. The losses are not specific to one cycle, and 
so any convenient cycle can be used for this study. For 
simplicity, the three-port flow divider cycle was chosen. In the 
flow divider, a single inlet flow is split into two outlet flows, one 
at higher stagnation pressure than the inlet flow, and the other 
at lower stagnation pressure than the inlet flow. No heat is 
added, so the apparatus is relatively simple (Fig. 1 ). 

In order to make an experimental study of losses, a wave 
rotor has been built at NASA Lewis Research Center, operating 
on the three-port flow divider cycle. This paper contains a brief 
statement of the philosophy of the experiment, a description of 
the experiment, and a summary of the measurements made. 
Finally, results are presented showing that reduction of the 
rotor-to-wall clearance gap leads to a large improvement in 
performance, and that friction and opening-time effects also 
play an important role in the performance. 


THE FLOW DIVIDER CYCLE 

The performance of a flow divider is conveniently indicated 
on a plot of the ratio of high stagnation pressure to inlet stag- 
nation pressure versus the ratio of low stagnation pressure to 
inlet stagnation pressure, with the mass flow ratio [3 as a 
parameter (Kentfield, 1 969). An upper limit to the performance 
of the flow divider can be calculated very simply, using what is 
called the acoustic approximation. In this approximation the 
following assumptions are made about the flow, 

(1) Flow conditions are constant within each region; regions 
are separated by waves. 

(2) Waves are not reflected at ports. 

(3) Waves travel at a single speed, which is the average of the 
wave propagation speeds on either side of the wave, and 
hence do not spread. 

(4) No change in entropy occurs across the waves. 

The results of this calculation are shown in Fig. 2(a). This 
performance, which is called the isentropic performance, is 
significantly higher than anything that can be achieved in 
practice, but it does illustrate the features of flow divider 
performance, i.e., that a large ratio of high pressure to inlet 
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pressure is only achieved at a low value of 3 - The dotted line 
terminating the curves of constant mass ratio 3 at the upper left, 
corresponds to an inlet Mach number of unity. In this approxi- 
mation, the curves of constant 3 are concave upwards, and 
maximum performance (i.e. largest high pressure ratio) will be 
achieved at an inlet Mach number of unity. 

A more accurate calculation can be made by using charac- 
teristics to evaluate the expansion out of the passages into the 
low pressure port, and including shock waves in the compres- 
sion portion of the calculation. This will be called the ideal cycle. 
The resulting flow divider performance is shown in Fig. 2(b). 
Now the curves of constant 3 are concave downwards, and the 
maximum performance (whether defined as pressure rise or 
efficiency) occurs for an inlet Mach number less than unity. A 
further, important, difference is where choking occurs. For the 
isentropic cycle, the inlet Mach number is greater than either of 
the outlet Mach numbers, and choking would occur at the inlet. 
In the ideal cycle, for values of 3 of 0.5 and less, reducing the 
low pressure leads to a Mach number of unity in the low pres- 
sure port while the inlet port is still subsonic. In other words, 
choking occurs at the low pressure port, not the inlet port. 
Although the ideal performance contains wave losses which are 
inherent in any real cycle, and so is a better approximation than 
the isentropic calculation, it can not be achieved in practise 
since instantaneous opening time and no friction are assumed; 
conditions which are obviously not attainable. 

The isentropic assumption is only valid for weak waves, 
which is the case for the lower right hand region of Fig. 2, 
where pressure ratios are close to unity. Comparison of 
Figs. 2(a) and (b) shows that the isentropic calculation works 
reasonably well for outlet Mach numbers M )0 less than about 

0.3. The experiments of Kentfield (1969) were entirely within 
this range, and his curves did not display a pressure ratio maxi- 
mum. For topping-cycle application, the pressure ratio should 
be as large as possible, and so one consideration in the pre- 
sent experiment was to operate the wave rotor under conditions 
for which the acoustic approximation would not be valid, and 
observe whether there is indeed a maximum in the perform- 
ance. Consequently, a design expansion ratio of 0.33 was 
chosen, corresponding to a Mach number of the expanded flow 
of 0.85. The expansion ratio is the ratio of the pressure in the 
low pressure port to the pressure in a rotor passage just before 
it is opened to the low pressure port. This ratio determines the 
Mach number of the flow in the low pressure port. For an 
expansion ratio of 0.33, the expansion wave exhibits significant 
spreading, as can be seen in the x-t diagram of the cycle shown 
in Fig. 3. 


LOSS MECHANISMS 

Although several wave rotors have been built in the 
past, only two studies of losses appear to have been reported, 
namely theoretical estimates of losses by Hoerler (1969) for the 
Comprex®, and by Kentfield (1963) for the flow divider. 
Kentfield gives the following losses as being in order of 
decreasing importance; 


1 . Basic wave effects 

2. Cell width, i.e. opening time effects 

3. Wall friction, passage entry and exit losses 

4. Leakage 

Hoerler calculated losses for the rotor itself, as percentage 
losses in efficiency for an experimental Comprex, and also for a 
hypothetical optimized Comprex. The results are given in 
Table 1. In addition, there are losses for ducts and stators 
external to the rotor. Hoerler states that leakage can dominate 
all other rotor losses, particularly for small machines, which is 
just the opposite of the conclusions by Kentfield. This disagree- 
ment reflects the different geometries of the two machines con- 
sidered, as will be explained below. 

It will be assumed that the friction, opening-time, and leak- 
age losses, will scale with one predominant parameter for each 
loss. The derivation of the parameter for each of these losses, 
plus the estimation of other losses will be outlined below. 

Basic Wave Effects 

Basic wave effects are the losses due to shock waves and 
spreading of expansion waves, which cause the difference in 
performance between the ideal cycle and the isentropic cycle. A 
real cycle can have shock waves instead of compression 
waves. There is a loss of stagnation pressure across a shock 
wave, leading to a reduction in performance. Expansion waves 
spread in space as they propagate, resulting in nonuniform 
velocities in the exit ports. Mixing of the nonuniform velocity dis- 
tribution to form a uniform distribution causes a stagnation pres- 
sure loss. These losses are unavoidable in a real device. The 
losses due to these effects depend on 3 and s, and are readily 
calculable if the timing is correct, i.e. if the ports open and close 
at the appropriate times as determined by the arrival, or launch- 
ing, of waves. For maximum performance, the timing will be dif- 
ferent for each value of 3 and s considered. An actual device 
will usually have timing fixed for one set of conditions, and 
operation off-design will create extra waves, causing additional 
losses. This can be handled computationally using a one- 
dimensional CFD code for wave rotor cycles (Paxson, 1995). 
The ideal cycle performance can not be measured directly in an 
experiment. It can be determined from experimental results by 
measuring performance as a function of the other losses, and 
extrapolating them to zero. 


Finite Passage Opening Time 

Since the passages have a finite width, there is a finite time 
taken for a passage to rotate past the leading or trailing edge of 
a port, and become fully open or closed. In the case of an inlet 
port, for which instantaneous opening of the passage would 
cause a shock wave to propagate down the passage, a finite 
opening time will result in a compression wave, which will 
steepen into a shock as it travels down the passage. The 
degree to which it steepens depends on the ratio of the opening 
time to the time taken by the wave to travel the length of the 
passage. Thus the relevant non-dimensional parameter is 
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1 - 0.2 


(4) 


Passage Opening Time ( JB \ ( al 

Wave Travel Time uj { l_j 


F = 0.32 — 
Dh 


LW 

v 


The speed of sound, a, will vary depending on which region of 
the cycle is being considered. For simplicity of definition, the 
inlet stagnation value will be used. Note that for any particular 
cycle, the rotor velocity U will be inversely proportional to the 
axial length of the passage L, and hence z is determined mainly 
by the value of the passage width, B. In addition to the effect on 
wave steepening, which may be advantageous, finite opening 
time will result in deleterious throttling losses when the passage 
is partially open. Further, finite opening time will create a grad- 
ual rise in velocity at the outlet ports, and finite closing time wilt 
create a gradual reduction in velocity at the outlet ports. The 
resulting non-uniform velocity distribution will result in a drop in 
stagnation pressure when it is mixed out to a uniform value 
downstream. Clearly this loss will increase as z increases. 


Friction 

Although the rotor passages are long and slender, their 
maximum length to width ratio is less than the entrance length 
for pipes. Thus the flow can be considered to have a boundary 
layer. This is confirmed by measurements of the radial velocity 
distribution in the high pressure port, showing a uniform velocity 
over the central 70% of the passage, and two-dimensional cal- 
culations of the flow in the entrance and high pressure ports of 
the experimental geometry (Welch and Chima, 1993) also 
showing a relatively small boundary layer. The flow in a 
passage open to the low pressure port can be thought of as 
flow over a flat plate, with the leading edge of the “plate" being 
the location of the leading edge of the expansion wave. A 
friction parameter can then be defined as the ratio of the drag 
force due to the flat plate to the product of inertial pressure and 
flow area, i.e. 

F = Z>(x) 1 pW 2 BH (2) 


The drag is related to the boundary layer momentum thickness 
at the end of the passage (Schlichting, 1979). Following Hoerler 
(1969), the boundary layer was assumed turbulent. With this 
assumption, the boundary layer momentum thickness is given 
by (Schlichting, 1 979) 

5 2 (x) = 0.036x(Wx/v)-°- 2 (3) 


The length of the boundary layer varies linearly with time as the 
expansion moves into the passage. By averaging over the time 
that the port is open, the friction parameter becomes 


This parameter will be taken as representative of the effect of 
friction on the cycle as a whole. The dominant factor in the 
friction parameter is L/D h . 


Heat Transfer 

Since it was not possible to measure heat transfer due to the 
lack of a sensor with a sufficiently rapid response rate, there 
seemed little point in making a separate assessment for heat 
transfer. It will depend on the same parameter as friction, so the 
two losses were simply incorporated together as one loss. 


Circumferential Velocity Loss 

If the angle a of the inlet duct is correct, the inlet flow will 
enter the passages smoothly at zero incidence in the relative 
frame (Fig. 4). This will only occur at design conditions (and 
even then, only if the design is correct!). Off-design, the flow will 
be at some angle of incidence, i, to the passage side walls, and 
work will be done on the gas, or extracted from it. This is 
equivalent to a flat plate compressor. The work done on the gas 
follows from the Euler equation for turbomachinery; 

C p AT = U(U-V in sina) (5) 


in which the stagnation temperature rise AT is defined as, 

AT = |3T Hj + (1 - P)T Lo - T jn (6) 


Measurements of motor power in some runs showed that the 
temperature rise could be as high as 20 °R, corresponding to a 
motor power of about 10 HP. Both the velocity and the flow 
angle a vary across the inlet port, and an overall value of the 
work done must be evaluated by integration. Depending on the 
angle of incidence, the relative stagnation pressure entering the 
rotor may be either greater or less than the absolute input stag- 
nation pressure. 

The wave action takes place entirely within the rotor, pro- 
ducing changes of pressure in the relative frame. However, the 
actual performance of a device is measured in the absolute 
frame, and will be affected by the rotor rotational speed. The 
experiments are run at three different rotational speeds, and so 
should be compared in the relative frame to remove the effect 
of the differing rotational speeds. Also, the ideal performance, 
as given above, is in the relative frame (since no correction for 
rotation was made, which is equivalent to assuming that the 
rotational velocity is zero). Thus it is necessary to correct the 
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observed absolute results to the relative frame. The flows leav- 
ing the rotor will have relative stagnation temperatures of mass - p^B (S W 0 W ^ + 5 p 0 p + 6 w 0 Woui )/w (11) 


Tr, = T, - U 2 /2C p (7) 


and hence the ratio of output relative to absolute stagnation 
pressure will be 


Defining a leakage parameter as the ratio of the mass leakage 
to the mass in a passage before it reaches the low pressure 
port, i.e. Ppag BHL, the leakage parameter becomes, if fa is 
assumed equal to p^. 


fs 

p, 


j _ 


Y/(Y-1) 


(8) 


In general, the relative pressure gain, Pr h /Pr ln , and the relative 
pressure loss, Pr^Pr,,,, will be greater than the absolute values, 
so that this effect is equivalent to a loss. However for the flow 
divider ( though not necessarily for other cycles), this effect is 
fairly small. 


G - 2(6 p 0 p + 5 w B w )/coHL (12) 


For operation of a specific wave rotor, wL is a constant. For 
convenience, relative values of 0 p , and 0 W , can be used such 
that 


0 P + ®w = 2 


(13) 


Leakage 

Leakage can take place radially from the passage to the 
casing if the passage is at high pressure, or from the casing to 
the passage if the passage is at low pressure. The result will be 
a "short-circuiting" from high pressure to low, leading to re- 
duced performance. In addition, circumferential leakage is pos- 
sible, from passage to passage. The pressure difference driving 
circumferential leakage is small except when a wave has 
reached the end of a passage. Thus circumferential leakage is 
likely to be small. For radial leakage, the rate of mass leakage 
will be proportional to the area available for leakage, which is 
26B at each end of a passage. The leakage will be into the 
passage while it is in a region of low pressure. From Fig. 3, the 
low pressure region on the inlet side extends over an angle 
0 , and on the outlet side, the low pressure region is the low 

pressure port, of extent 0 P , and the region between the low 
pressure port and the high pressure port, of extent© .Thus 

w out 

the time spent at low pressure is 


(9) 


at the inlet side of the wave rotor, and 


'■V“* e wJ“ do) 


at the outlet side of the wave rotor. The amount of mass leaking 
into a passage will be 


and the proportionality constant chosen so that the leakage 
parameter is defined as; 

G = (6 p 0 p + 6 W 0 W )/H (14) 

Then if the leakage gaps are equal, i.e. 5 P = S w = 6, this re- 
duces to 

G = 25/H (15) 


Based on this derivation, one would expect that the leakage in 
the Comprex studied by Hoerler (1969), would be similar to that 
in the experiments of Kentfield (1963), since the ratio 25/H was 
virtually the same for both devices, but Hoerler claimed leakage 
was important, whereas Kentfield said that it was not. However, 
the two devices operated on quite different cycles, and whereas 
the experiment of Kentfield used a similar cycle to that of the 
present experiment, the Comprex cycle was different. In the 
Comprex, which used a four-port cycle, the inlet port and 
exhaust port were at approximately the same pressure, which 
was significantly lower than the two high pressure ports. The 
cavity was therefore likely to be at some intermediate pressure. 
This means that the assumption above that is prob- 

ably not correct for the Comprex, and also there would be leak- 
age into the inlet port. Both of these effects would lead to more 
leakage. An additional difference was the relative extent of the 
low pressure regions, with more of it in a port region for the 
Comprex. It thus seemed important to ascertain in the present 
experiment whether leakage at a wall is more, or less, important 
than leakage at a port. In order to do this, the wall gap and port 
gap were varied independently. Thus separate port and wall 
leakage parameters were defined, i.e. 
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G p = W H 

(16) 

G„ = <WH 

(17) 


Relative values of 6 P = 0.71 and 0 W * 1 .29 were used, so that 
when the wall and port gaps are equal, the total leakage 
parameter, 

G = G p + G w =(0.715 + 1.295)/H = 25/H (18) 


in agreement with Eq. (15) above. 


Entrance and Exit Losses 

The wall between the passages has a finite thickness. There 
is consequently an area change between flow in the ducts, and 
flow in the passages. The squared ends of the walls will present 
an obstruction to the flow on entering, and a drag on leaving, 
leading to stagnation pressure losses. This effect was calcul- 
ated by assuming the leading edge is at the entering stagnation 
pressure, and the trailing edge is at the downstream static 
pressure. The resulting force on the flow is included in the mo- 
mentum equation, which is then solved with the continuity and 
energy equations to give the downstream stagnation pressure. 
The losses in stagnation pressure are less than 2% for the 
narrow passages, and less than 1% for the wide passages. 


EXPERIMENTAL DESIGN 

In order to be able to obtain experimental values of the 
losses due to opening time, friction and port and wall leakage, 
the parameters x, F, G p and G* must be varied. An efficient way 
to formulate an experiment to obtain empirical fits to data when 
there are three or more variable parameters is the Box- 
Behnken scheme (Box and Behnken, 1960). This scheme is 
illustrated for three variables in Fig. 5. Imagine a box, each side 
of which extends from the minimum value of the corresponding 
variable to the maximum value of the variable. Experimental 
readings are taken at the points indicated in the middle of each 
side, together with three replicate points at the center of the 
box. It is then possible to fit the results with a second-degree 
polynomial. For example, if the measured variable is the effi- 
ciency n, and the independent variables are t, F, and G, then 
the fit will be of the form 

H = b 0 + b,x + b 2 F + b 3 G + b u x^ + b^F 2 + b^G 2 + b 12 xF 

+ b 13 xG + b^FG (1 9) 

where the constants b 0 ,b„ and b„ are determined from the ex- 
perimental measurements. The replication of the center point 


provides an estimate of the experimental error. The scheme for 
four variables is equivalent, but harder to illustrate since draw- 
ing a four-dimensional box is rather difficult! The set of runs that 
must be made in order to evaluate the constants was deter- 
mined using commercial software (Seshadri and Deming, 
1990), and is listed in Table 2. The same software was also 
used to determine the constants from the experimental data. 

The present wave rotor experiment was formulated as a four 
variable design, with t, F, G p , and G* as the independent vari- 
ables. The dependent variable, or response, was taken to be 
the efficiency, which has been defined for a flow divider 
(Kentfield, 1969) as 


(3 

jv 

(1 - 3) 



In designing the experiment, the actual rotor dimensions had to 
be determined. A review of the literature showed that the geom- 
etry had differed significantly for rotors built in the past, some- 
times with no indication of the value of the end-wall to rotor 
gaps. The details of several past rotors are given in Table 3. 
The non-dimensional opening time t has varied over a range of 
0.1 to 0.35. The major component of the friction parameter, 
namely the ratio L/D h , has varied between 3.7 and 84. Finally, 
the leakage parameter G (previous workers have not differ- 
entiated between G p and G*,), has varied from 0.005 to 0.017. 
The objective in designing this experiment was to cover as 
much of this range as possible. In order to do this, the actual 
experiment consisted of two different rotors, both 12" in diam- 
eter, but one 18" long, the other 9" long. Both were built with 
passages 0.25" wide, and 0.4“ high, with 120 passages per 
rotor. After a series of runs at 0.25" passage width, every other 
wall was removed, and another series of runs was made at 
approximately twice the passage width. Ducts were designed 
for operation at n = 1.5 for both the 9" and 18" long rotors. 
Adding a different set of ducts for operation at n = 3 with the 1 8“ 
rotor, provided a combination giving three values of t, and three 
values of L/D h , covering the range 20 to 58, and hence three 
values of F. The value of n = 1 .5, rather than n = 1 , was used to 
provide a long region between closing of the high pressure port, 
and opening of the low pressure port, to allow any remaining 
waves to die out, and give the uniformity prior to opening the 
low pressure port that was assumed in calculating the cycle. In 
fact, only one actual cycle was used; what is meant by n = 1 .5, 
and n = 3, is the number of cycles there could be per revolution 
at the timing used for the one cycle. In order to vary the port 
leakage gaps, the ports were built as inserts, supported on a 
flange. Placing shims under the flange permitted variation of the 
rotor to port gaps. The rotor to wall gap was varied by providing 
moveable end-walls at each end, consisting of disks, which 
could move axially. Springs at three locations pushed the disks 
away from the rotor, and three screws forced the disks towards 
the rotor. Adjustment of the screws gave the desired gap spac- 
ing. The minimum spacing that could be used safely was 
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0.005", larger than had been hoped for. This gave a value of 
G = 0.025, which was greater than the values used by other 
workers due to the small passage height of the present 
experiment. 

The experimental results were measured in the absolute 
frame of reference. Since in the experiment, different rotational 
speeds, and different ratios of wall thickness to passage width 
were used, the circumferential loss and entrance and exit 
losses were not the same from run to run. Consequently the 
absolute pressure ratios were converted to the relative frame, 
and corrected for entrance and exit losses. The differences 
between the corrected pressure ratios as the geometry is varied 
should then be the differences due to the variation in opening 
time, frictional, and leakage losses. 

A photograph of the apparatus, with the 18" rotor in place, 
is shown in Fig. 6. 


EXPERIMENTAL MEASUREMENTS 

The efficiency of the flow divider (Eq. 20) is determined by 
the ratios of P h /P in and P to /P in at a particular value of (3. Higher 
values of both ratios correspond to higher efficiency. The nec- 
essary measurements are therefore the mass flows in each 
port, and the stagnation pressures of the inlet, high pressure, 
and low pressure flows. The mass flows were measured with 
standard orifice meters. The ports are obviously sections of an 
annulus. A transition piece in front of the inlet port took the flow 
in the inlet pipe, and converted it to the port shape, accelerating 
the flow in the process, as well as bringing it onto the rotor at 
the correct angle. Immediately upstream of the transition piece, 
but downstream from the orifice, was a diagnostic spool, with 
three wall static taps, five pitot tubes, and a thermocouple. The 
inlet stagnation pressure and temperature were determined by 
measurements at this spool. Similarly, the output ports had 
transition pieces to take the flow from the port shape back to 
round, and which also acted as diffusers. The downstream area 
of these diffusers was fixed by the exhaust pipe diameter. With 
the diffuser area ratio determined, the length was chosen to 
give maximum diffuser efficiency, using the diffuser perform- 
ance curves of Mattingly et al. (1987). A diagnostic spool was 
placed immediately at the exit of each diffuser, with the inten- 
tion of using the measurements to evaluate exit stagnation 
pressures. However, the velocity distribution was found to be 
very non-uniform, and it was not clear whether the resulting 
stagnation pressure would be reliable. In addition, measure- 
ments made at the ports were used to calculate a stagnation 
pressure. There were 5 static pressure taps on the top and 
bottom of each port, and four pitot tube installations. The pitot 
installations carried either a rake of 5 pitot tubes to determine 
radial velocity distribution, or a tube-type combination probe 
(Glawe and Krause, 1974), to determine center-line velocity and 
direction. In the high pressure port, velocities are low, and 
relatively uniform, and an average of the individual port stagna- 
tion pressure measurements was taken as the port stagnation 
pressure. In the low pressure port, the velocity distribution is 
very non-uniform, both radially, and circumferentially. The 
measurements were used to create a circumferential velocity 


distribution, and also, by using a rake probe, to create a radial 
velocity distribution. These distributions were used to create 
mass, momentum, and energy integrals for a mixing calculation 
(Foa, 1960), from which the stagnation pressure for a uniform 
downstream velocity was evaluated, and used as the low pres- 
sure port stagnation pressure. The pressure measurements 
were steady state. For some runs, a dynamic pressure trans- 
ducer was installed 4" downstream of the low pressure port. 
The signal from this transducer was a sine wave at the passage 
passing frequency, with a peak to peak value of 0.5 percent of 
the steady state pressure. Thus the pressures were essentially 
steady state. 

All steady state pressure measurements were recorded 
through an electronically scanned pressure (ESP) measure- 
ment system. The pressure measuring system automatically 
self-calibrates every 20 minutes to maintain a 0.1% accuracy. 
The ESP system communicates through an Institute of Electri- 
cal and Electronics Engineers 488 interface to a state of the art, 
real time data acquisition system designed at NASA Lewis 
(Fronek et al., 1987). The data system was designed for small 
to medium sized aeronautics facilities, most of which are cur- 
rently testing rotating machinery. The system acquires data, 
converts it to engineering units, computes test dependant per- 
formance calculations, and displays the information in alpha- 
numeric or graphical form. The cycle time is one second. 

Control of the wave rotor flows was by a butterfly valve in 
each leg. The supply pressure was around 55 psia. The inlet 
valve was adjusted to throttle this supply pressure down to the 
desired inlet stagnation pressure, usually 30 psia. The expan- 
sion ratio was set by the low pressure valve, based on a low 
pressure port static pressure reading, and a pressure tap in the 
endwall giving the passage pressure just before opening to the 
low pressure port. The high pressure valve controlled the mass 
ratio 3» based on the input and high pressure orifice mass flow 
readings. The rotor was turned by a variable speed electric 
motor with a constant speed control. An independent measure- 
ment of the rotor speed was also made. 


RESULTS 

For the statistical experiment, runs were made at a 
nominally constant value of 3 = 0.37, mostly at an input stagn- 
ation pressure of 30 psia, varying e from 0.33 to 0.8, unless the 
full range was not accessible. Because it was not possible to 
achieve 3 = 0.37 exactly on every run, runs were also made at 
a nominal 3 = 0.36, and 3 = 0.38. From these extra runs, a local 
value of the derivatives of P h /P in an F} 0 /P in with respect to 3 
could be calculated, and used to correct the runs at nominal 
3 = 0.37 to values corresponding to exactly 3 = 0.37. In 
addition, corrections for conversion to the relative frame, and 
for entrance and exit losses were made. Each run was dupli- 
cated, and then the average of the corrected values of P hi /P jn 
and Pn/Pj,, for the two runs was taken as the final result. 

The results were plotted as P h /P in versus PJP^. A sample of 
the results showing performance changes for the different 
rotors with both port and wall gaps set equal to 0.010", is given 
in Fig. 7. All these runs were made at an inlet stagnation pres- 
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sure of 30 psia. The runs made with the 18" rotor at n = 3 have 
a larger opening time than the runs with n = 1 .5, but the same 
value of friction factor. Thus it is seen that increasing the open- 
ing time reduces performance. Runs with the wider passages 
are superior in performance to runs made with the narrow pass- 
ages. The runs with the wider passages have a larger opening 
time, but lower friction factor than the runs with the narrow 
passages. For these runs then, reduction in friction is more 
beneficial than the increase in opening time is deleterious. This 
is more evident for the runs with the 18" rotor than it is for the 
runs with the 9" rotor, indicating that friction is becoming less 
important for the shorter rotor. 

Results obtained with the 18" rotor at three different gap 
spacings are shown in Fig. 8. The three runs at a gap spacing 
of 0.010" are the replicated runs at the center point, giving an 
idea of the experimental error due to lack of reproducibility of 
gap setting. It is clearly much less than the change in perform- 
ance caused by a gap change of 0.005", so the observed effect 
is real. The variation at a fixed gap setting is of the order of 
0.003, much less than the error due to gap setting, showing that 
irreproducibility of gap setting is the major source of experimen- 
tal error. There is a significant increase in performance as the 
gap spacing is reduced. 


EVALUATION OF THE RESULTS 

At each value of expansion ratio, the corrected values of 
Ph/Pin and P to /F? n were used to calculate an efficiency, and this 
was entered into the statistical program. The program then 
evaluated the values of b 0 , b, and b H , and provided the confid- 
ence estimate of each value. Values with low confidence were 
dropped from the model. For example, at an expansion ratio of 
0.6, (which is intermediate between the maximum pressure 
ratio at an expansion ratio of 0.55, and the maximum efficiency 
at an expansion ratio of 0.65), the final model was 

H = b 0 + b 3 G p + b 4 G w + b^x 2 + b^F 2 + b 12 xF (21 ) 

The values of the constants are given in Table 4. This remark- 
ably simple model shows that the efficiency decreases linearly 
with the end-wall to rotor gap spacing. Moreover, since lb 3 l is 
larger than lb 4 l, when the relations in Eqs. (16) and 17) are 
entered for the leakage parameters, the efficiency is almost 
equally sensitive to port gap and wall gap, 

H = 0.42 - 0.845/H - 0.725,/H - 1 .3t 2 - 1 .74F 2 + 0.70xF (22) 

The linear relationship between efficiency and leakage param- 
eter is confirmed in Fig. 9, showing efficiency versus leakage 
parameter for the 9" rotor. This result indicates that values of 
leakage parameter of 0.005 or less are desireable for good 
efficiency. 

The model can be used to extrapolate the results to zero 
leakage parameter. A contour plot of efficiency at no leakage, 
as a function of opening time and friction factor, derived with 
the full model, is given in Fig. 10. It will be seen that there is a 
ridge of high efficiency sloping away from a maximum at the 


origin, and falling off rapidly along the axes. This result can be 
used to derive an optimum passage width for a given wave 
rotor. Assume that the rotor length and passage height are 
fixed, and that the ratio of B/H is 2 at a friction factor of 0.4. As 
the value of B is decreased, the opening time gets less, and the 
friction factor increases. As B is increased, the opposite is true. 
Thus the locus of points along a curve of constant rotor length, 
but varying passage width is approximately a hyperbola. The * 
exact curve is plotted in Fig. 10, labelled "L = constant". This 
curve is tangent to the contour of efficiency = 0.48 at the point 
F = 0.42, x = 0.185. As will be seen, this is also the point of 
highest efficiency for this curve, and therefore the passage 
width appropriate to this point is the optimum for this rotor. 
Examination of the other contour plots indicates that the 
optimum for different length rotors will lie along the line 

F = 2.25x (23) 

The contour plot shows that the maximum efficiency is at the 
origin, i.e. as small a value of friction factor and opening time as 
possible. This may not be the case in reality, for two reasons. 
First, the model is strictly valid only in the space covered by the 
experiment, which was the region x = 0.085 to 0.37, F = 0.37 to 
1.12. The origin is outside this region, and extrapolation of the 
model to the origin may give erroneous results. Secondly, the 
results have been referred to the relative frame. As one gets 
closer to the origin, the rotor length gets smaller, and the rotor 
speed goes up. This can result in the absolute performance 
being significantly below the relative performance, so that a 
longer, slower, rotor may have better performance. Such an 
effect was found for four-port wave rotors in an optimization 
study (Wilson and Paxson, 1995). If the extrapolation to the 
origin is correct, the result should be the efficiency calculated 
from the ideal cycle, or rather, since an expansion ratio of 0.6 is 
off-design, the efficiency calculated with the one-dimensional 
model of Paxson (1995). Extrapolation with the full model gave 
an efficiency of 0.54, compared with an efficiency of 0.56 
calculated with the one-dimensional model. 


CONCLUSIONS 

Operation of a three port wave rotor has shown that, at a 
constant value of 3 = 0.37, the maximum pressure ratio occurs 
at an expansion ratio 0.55, and the maximum efficiency at an 
expansion ratio of 0.65. These expansion ratios are significantly 
higher than the lowest possible expansion ratio. The loss of effi- 
ciency in a wave rotor due to leakage has been demonstrated 
experimentally to be linear in the end-wall to rotor gap spacing. 
Values of the leakage parameter G = 25/H less than 0.005 are 
required for good efficiency. Friction and opening time losses 
cause the efficiency to decrease with the square of the friction 
factor, and the square of the opening time, with an offsetting 
term increasing as the product of the friction factor and the 
opening time. Optimum wave rotor efficiency will be obtained if 
the rotor is designed so that the friction factor is 2.25 times the 
dimensionless opening time. 
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TABLE 1.— LOSS OF EFFICIENCY IN A COMPREX® 



Source of loss 

Loss in experimental 
COMPREX 

Loss in optimized 
COMPREX 

Shocks and Fans 

4.5% 

3.3% 

Passage opening/closing 

4.1% 

4.2% 

Friction (turbulent) 

7.2% 

2.1% 

Heat transfer 

3.3% 

2.1% 

Rotational velocity 

4.2% 

2.0% 

Leakage 

10.1% 

1.1% 

Mechanical 

2.3% 

0.1% 


interface mixing 


Total 



TABLE 2.— SET OF RUNS FOR THE STATISTICAL EXPERIMENT 


Run Rotor Passage 

number length, width. 


in. in. 



Note that runs 1 1 , 1 6, 23 are replicates. 
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TABLE 3. — COMPARISON OF DIFFERENT WAVE ROTORS 



Thayer 

Weatherston 

Pearson 

Hoerler 

Kentfield 

This work 

Machine type 

Pressure 

exchanger 

Pressure 

exchanger 

Wave 

turbine 

Pressure 

exchanger 

Flow 

divider 

Flow divider 

Length- L' 

15.8 

66 

3.5 

4.25 

11 

9.18 

Diameter-D" 

14.4 

60 

9 

3.23 

8 

12 

Passage width-B' 

0.4 

0.55 

0.7 

0.19 

0.66 

0.25, 0.54 

Passage height-H' 

1.5 

1.43 

1/5 

1.1 

2.2 

0.4 

Hydraulic dia-D" h 

0.63 

0.79 

0.95 

0.32 

1.05 

0.31,0.46 

RPM 

1 960 

2 700 

18 000 

11 200 

5.500 

1850-7400 

Cycles/revolution n 

2.5 

1 

1 

2 

3 

1.5,3 

Leakage gap-5 

0.004-0.013 



0.004 

0.007 

0.005, 0.01,0.015 

Ratio UD 

25 

84 

3.7 

13 

10.5 

20-58 

Ratio 25/H 

0.005-0.017 



0.007 

0.006 

0.025-0.075 

X 

0.2 

0.1 

0.3 

0.33 

0.35 

0.08-0.35 


TABLE 4. — VALUES OF THE CONSTANTS IN EQUATION (13) 


Constant 

EH 

m 

b 4 

b„ 


bi2 

Value 

0.42 

-1.18 

-0.56 

-1.30 

-1.74 

0.70 




































































































Upstream 
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Pitot probes-either 
0 rakes or directional 
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pressure 
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mass/total 
head 
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Pressure 
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in endwall 


Dynamic pressure 
transducer - on rotor, 
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Figure 1 . — Schematic diagram of the apparatus with the 
diagnostics used. 


High pressure 
X 



Low pressure 
outlet i 


Main expansion wave 

Figure 3. — Characteristics calculation of the flow divider 
cycle for € = 0.33. 



P|</^in P|</Pjn 


Figure 2. — Performance of a three-port flow divider calculated for (a) the isentropic 
cycle (left), and (b) the ideal cycle. 


12 






Figure 5. — Illustration of the Box-Behnken design of an 
experiment for three variables. Runs are made at values 
of the variables corresponding to the solid circles shown. 
In addition the center point, indicated by the three 
concentric circles is repeated three times. 



Figure 6. — A photograph of the three-port flow divider apparatus. 
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03 0.4 05 0.6 0.7 03 


Pio/Pin 

Figure 7. — A summary of the experimental results, showing 
the variation in performance with changes in opening time 
and passage length. 



Leakage parameter G 

Figure 9. — Experimentally observed flow divider efficiency 
versus leakage parameter for the 9" rotor, at an expansion 
ratio of 0.6. 


P h/Pin 


Friction factor, F 



P|c/Pin 


Figure 8.— A summary of the experimental results obtained 
while varying the leakage gap. Results for a gap of 0.01 0" 
represent the center point, which was run three times. 



o 0.1 0.2 03 0.4 
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Figure 1 0. — Contour plots of flow divider efficiency as a 
function of opening time and friction factor, calculated 
with equation 22. 
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